function z = rnd_mvn(mu,sigma)
z = mu + chol(sigma)'*randn(length(mu),1);
